function thetaFunc = thetaFunc(vec1,vec2)
    
    thetaFunc = acosd(dot(vec1,vec2)/ ... 
        (norm(vec1)*norm(vec2)));

end